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Massive gravity, galileon and braneworld models that modify gravity to explain cosmic acceler- 
ation utilize the nonlinear field interactions of the Vainshtein mechanism to screen fifth forces in 
high density regimes. These source-dependent interactions cause apparent equivalence principle vi- 
olations. In the weakly-screened regime violations can be especially prominent since the fifth forces 
are at near full strength. Since they can also be calculated perturbatively, we derive analytic so- 
lutions for illustrative cases: the motion of massive objects in compensated shells and voids and 
infall toward halos that are spherically symmetric. Using numerical techniques we show that these 
solutions are valid until the characteristic scale becomes comparable to the Vainshtein radius. We 
find a relative acceleration of more massive objects toward the center of a void and a reduction of 
the infall acceleration that increases with the mass ratio of the halos which can in principle be used 
to test the Vainshtein screening mechanism. 



I. INTRODUCTION 

In models that seek to explain the current acceleration 
of the cosmic expansion by modifications to gravity or the 
addition of universal fifth forces on cosmological scales, 
such modifications must be hidden from local tests by 
so-called screening mechanisms. Screening mechanisms 
invoke nonlinearity in the equations of motion for the 
field that mediates the extra force. 

The Vainshtein screening mechanism p] [2] was first in- 
troduced in the context of massive gravity to suppress the 
propagation of additional helicity modes [3J 0]. Here non- 
linear derivative interactions of the field act to screen the 
fifth force within the so-called Vainshtein radius around 
a matter source. The Vainshtein mechanism occurs not 
only in modern incarnations of Boulware-Deser [5j ghost- 
free massive gravity [5HS] but also in galileon cosmology 
PMil] and braneworld models [T7H23]. 

In these models, all bodies accelerate equivalently in 
the total field of the fifth force and hence obey a micro- 
scopic equivalence principle. Nonetheless screened bodies 
do not move as test bodies in an external field leading to 
an apparent or macroscopic violation of the equivalence 
principle |24j . Instead, nonlinearity in the field interac- 
tions causes interference between the external and body 
field in forming the total field [25] . In the Vainshtein 
mechanism, this interference occurs when the second spa- 
tial derivatives of the fields becomes large enough that 
the self-interaction terms become important, i.e. within 
the Vainshtein radius of the external and body sources. 

A general technique was recently introduced for deter- 
mining such effects in two-body systems by considering 
the effective density generated by the nonlinear inter- 
action [26 . It was applied to the Earth- Moon system, 
which exhibits strong screening since the orbit of the 
Moon is well within the Vainshtein radius of both the 
Earth and the Moon. Violations of the equivalence prin- 
ciple thus appear as a small mass-dependent correction 
on the already suppressed effect of anomalous perihe- 



lion precession. Here we study the weakly-screened limit 
where the Vainshtein mechanism is only beginning to op- 
erate. These cases have the advantage that fifth forces 
are only weakly suppressed and that equivalence princi- 
ple violations are analytically solvable by a perturbative 
expansion. Weak screening is applicable to cosmological 
situations such as voids and the outskirts of dark matter 
halos. 

The outline of this paper is as follows. We begin in fjll] 
with a review of the Vainshtein mechanism and develop 
the effective density approach for the weak-screening 
regime. In §III[ we consider several examples where the 
acceleration of massive bodies can be analytically cal- 
culated. In j jIV| we test the limits of validity for the 
weak-screening approximation. We discuss these results 
in|Vl 

II. WEAK SCREENING 

For models that exhibit Vainshtein screening, the or- 
dinary Newtonian potential is modified by the addition 
of a scalar field <f> which itself obeys a nonlinear Poisson 
equation 

3/3V 2 + N[(j>, <)>] = 8nG5p, (1) 

where ft is a parameter that determines the cou- 
pling to matter density fluctuations from the cosmo- 
logical mean Sp = p — p. For definiteness, we take 
the Dvali-Gabadadze-Porrati (DGP) braneworld exam- 
ple |27j where the derivative operator is given by 

N[<f>!,<h] = r\ [V 2 ^V 2 2 - ViV^iVV> 2 ] (2) 

in the quasistatic approximation. Here r c is the crossover 
scale and determines the strength of the nonlinear inter- 
actions. We have written out this derivative interaction 
term in a general bilinear form to facilitate its use in 
two-body systems [2"S|. 



For a single, spherically symmetric body of mass M it 
is straightforward to show that fifth forces are screened 
within the Vainshtein scale (e.g. [2"5] ) 
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For example, we can define a weak-screening regime ex- 
terior to this scale where 126 
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For two bodies, there is an additional screening that 
can be expressed as the interference between the fields of 
the individual bodies. Given two sources SpM and Sp m , 
which individually produce fields 4>m and 4>m satisfying 



3f3V 2 4> M ,m + N[(p M>m , 4>M,m] = ^G8p M ,r 
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jointly <j) 

m ~t~ no longer solves Eq. ( JT| ) for Sp^i -\- &Pm- 
Instead, the interference field 126 



solves Eq. ([I]) if 
3/3V 2 0a = 



-2N[<j> M ,<f>r, 

-N[<t> A ,<p A ] 



2N[4> M + 4>m,4>l 



(6) 



(7) 



Thus the real density fields 5pm and Sp m are eliminated 
in favor of an effective density field given by the nonlinear 
interaction terms. Note that the real density fields can 
themselves be composite systems of n individual bodies 
so long as the solutions 4>m and <j> m are known. In gen- 
eral this system remains a nonlinear Poisson equation, 
amenable only to numerical techniques [3TJ [551 HI] • 
Now let us consider the weak-screening limit. By tak- 
0, we can iteratively solve for joint screening 
effects order- by-order in r 2 . To leading order </>a = 0{r 2 ) 
and Eq. Q can be approximated as 
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which is now a linear Poisson equation with an external 
effective density. Using this leading order expression, we 
can find the first order correction by solving 
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etc. Convergence of this series checks the validity of 
the approximation. By dimensional analysis, the series 
should converge so long as t*m/s "C 1 and r* m /s <C 1 
where s is the typical physical scale of the system, e.g. 
the separation between bodies (see {IV). 

The weak-screening limit for two-body interactions can 
apply even if near the individual bodies the individual 
fields enter a strong-screening regime. For example near 
the location of mass m, <f> rn may actually be strongly 
self-screened if the physical size of the body is smaller 
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FIG. 1. Weak screening examples: (a) compensated shells 
with 8pM — in the interior and mass ±M in the shells; (b) 
compensated void with an underdensity of mass — M in the 
center and +M in the shell; (c) infall toward onto a spherical 
mass M. The smaller point particle, whose motion we study, 
is labeled m. 



than its Vainshtein radius. All that is required for two- 
body weak screening is that the interference source N is 
dominated by the weakly-screened far-fields and that it 
generates an interference field (f>A that is nearly a pure 
gradient across the Vainshtein radius of m [23]. Higher 
order contributions from Eq. ^ are therefore small since 
the second derivatives of </>a are small where those of <j> m 
are large. 



III. ANALYTIC EXAMPLES 

In this section, we consider three examples of equiv- 
alence principle violation in weakly-screened systems 
where the mass-dependent acceleration of bodies in the 
system can be solved analytically: empty, compensated 
mass shells; negative density fluctuation, compensated 
voids; and infall into a spherically symmetric body. 



A. Shells 

In Newtonian mechanics the shell theorem says that 
a body inside a spherically symmetric shell does not ex- 
perience a force from the shell. Despite the lack of a 
(9) purely 1/r 2 force under the Vainshtein mechanism (see 
e.g. Eq. [4]), it remains true that a test body of infinitesi- 
mal mass inside the shell does not experience a force (see 
e.g. [23] ) . However for a finite mass m, the nonlinear in- 
teraction between the body field and shell field causes a 
force unless the body is at the center of the shell. This 
effect is common to gravitational field equations that are 
nonlinear. In general relativity its impact is suppressed 
by the strength of the nonlinearity, i.e. the ratio of the 
Schwarzschild radius of the body to the separation Gm/s, 
and in Modified Newtonian Dynamics it appears unsup- 
pressed in the deeply nonlinear regime [30] . 
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FIG. 2. Effective density as a function of z = rcosO for the 
shells case of Fig.[I^ with = 1.2r\ and — > oo and several 
choices of the displacement d of the particle from the center of 
the shell (crosses). As the particle is displaced further along 
the +z axis, the effective density in the shells exhibits an 
increasing asymmetry which leads to the net force. 



To set up this first test case, we take bpu to contain 
a constant-density shell that encloses a mass of +M be- 
tween ri < r < r- 2 - In order to ensure compatibility 
with cosmological boundary conditions where the mean 
density fluctuation is zero, we take an additional outer 
shell of — M between r 2 < r < r 3 . Note that in the limit 
r 3 — > oo the result will be the same as an uncompen- 
sated single shell. We obtain the 4>m solution to zeroth 
order in r 2 by superimposing the Newtonian solutions for 
a tophat constant density enclosing a mass M t h within 
the radius 
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Superimposing tophats of the appropriate mass, 4>m = 
const, for r < ri and for r > r 3 . 

Next we take a small body of mass m displaced along 
the z axis at z — d < r\ so as to be inside the shell 
(see Fig. []J,) . It suffices to consider this mass to be un- 
compensated. If we place a compensating shell of — m 
at r > r 3 , we obtain identical results since there can 
be no interference with the cf>M — field there. More 
specifically, N[<j)M,<j>m\ has compact support and is only 
non-vanishing for r\ < r < r 3 



N[<f> M ,<f> m ] = cj f rlF(r,e) 
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where r* m is the Vainshtein radius of m in isolation. 

Even if cj> m were in the deep screening regime near the 
location of to, there would still be no effective density 
there since 4>m is constant inside the shell. Thus we can 
safely take the limit of a point mass particle. In Fig. [2] 
we show how the effective density changes as the point 
mass is taken from the center of the shell toward the edge. 
Note the asymmetry that develops between the near and 
far side of the shell. 

We can build our solution for 4>a from that of the fol- 
lowing electrostatics-like Poisson equation 



V 2 $ = 
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r < r < oo 
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More specifically, we are interested in calculating the po- 
tential gradient V$ = d$/dz at a position of the particle 
r = d and 9 = 0, or the "electric field" at that point 
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In the d — > limit, x — > oo and I(x) 

We again superimpose these solutions with the appro- 
priate tq to construct the full solution 



dl(ri/d)-l(r 2 /d)} 
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Eugeny Babichev Note that c g scales as the shell acceler- 
ation GM/f3rf suppressed by the cube of the Vainshtein 
scale over the characteristic distance. 

Examples of the acceleration as a function of djr\ are 
shown in Fig. [3j One interesting limit is an infinites- 
imally thin shell r 2 — > r\ with vanishing impact from 
cosmological compensation r 3 — > oo 
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Consider the cases where the particle is near the center or 
near the shell. For the former case, there is a suppression 
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FIG. 3. Acceleration or field gradient as a function of distance 
d/ri from center for various shell widths given by r%jr\ and 
r% — ¥ oo. As the particle approaches the shell, the acceleration 
increases as s -3 , where s = ri — d, with a proportionality 
constant that reaches a maximum as the shell width goes to 
zero r%jr\ — > 1. 



of ( r *m/ r i) 3 (d/ri) from the scale of the Newtonian ac- 
celeration of the shell mass. The first factor accounts for 
the weak nonlinearity of the system. The second factor 
accounts for the fact that by symmetry the force must 
vanish if d = 0. 

In the opposite regime where the particle is closer to 
the shell than the center, the acceleration scales strongly 
with distance to shell s = r\ — d. In particular as the 
particle approaches the shell 
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such that there is a suppression factor of the Vainshtein 
radius compared to the distance to the shell (see Fig. [3]) . 

In this case, we see the impact of the nonlinearity of 
the field equations unobscured by symmetry. Once the 
Vainshtein screened regime of m starts to overlap the 
shell, the field of the shell effectively creates a boundary 
condition there, much like the establishment of a mirror 
charge by a conductor. Since the Vainshtein radius of 
the particle depends on its mass, the equivalence prin- 
ciple is macroscopically broken: for an attractive fifth 
force, the acceleration of a particle towards the center 
will grow linearly with the mass m. The sign of the effect 
is determined by the fact that the Vainshtein correction 
in Eq. Q makes the force fall away with distance less 
rapidly than 1/r 2 and so the far side of the shell that 
contributes more force than the near side. 



B. Void 

A closely related and more observationally relevant 
case is a cosmological void. In this case instead of an 
interior at the mean cosmological density and 5pM = 0, 
we have an underdensity. We can idealize the void 
as a spherical tophat of spatially constant underdensity 
&Pm — pAy where — 1 < Ay < for r < r% with a total 
negative mass fluctuation of 
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The void is surrounded by a positive density shell of mass 
+M for r\ < r <r2 such that the total mass fluctuation 
exterior to T2 is zero (see Fig. JTb). The Vainshtein scale 
of such a void is 
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and so for H$r c / ' {3 < 1, a void of any size r\ is in the 
weak-screening regime throughout its interior. 

Despite the fact that the interior now has a finite SpM, 
the void and shell case are nearly identical. The reason 
is that if 4>m is the r 2 -field associated with a constant 
density source, then 



N[(j> M ,(l>m} = ^V^mV 2 ^ 
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regardless of the form of 4> m [2"fJ] . 

For a particle of mass m, the effective density therefore 
vanishes inside the void except at the position of the test 
mass. If we take the particle to be a constant density 
spherical tophat with p m that vanishes outside its radius 
r m < n 



N[4>M, <t>m] = 
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The r < r± term provides no acceleration at the po- 
sition of the test mass by symmetry and the shell term 
gives 
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which has a form very similar to Eq. (17 1. In particular 
the infinitesimal width r 2 — > r\ case has the limiting 
behaviors 
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where recall s = r\ — d is the separation between the 
particle and the shell. In the later limit, the shell and 
void cases are identical as one might expect. 



4 



Halos in compensated spherical voids therefore accel- 
erate differently depending on their mass when they are 
near a Vainshtein radius r* ro of the shell. The Vainshtcin 
radius of a halo scales with its virial radius. Taking the 
virial radius of the halo r # as the radius out to which the 
average interior density reaches Ah ~ 200, we obtain the 
Vainshtein radius of a halo of mass m 
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For Hor c /f3 < 1, halos that are near a virial radius of 
the compensating shell will be accelerated away from the 
shell. 



C. Infall 

Lastly we consider the problem of two point masses, or 
tophat constant density bodies of radius negligible com- 
pared with their separation. This problem was solved 
numerically in the strong-screening regime in Ref. |26j . 
Here we consider the weak-screening regime as might be 
appropriate for a satellite dark matter halo of mass m 
falling through the virial radius of a parent halo of mass 
M. 

In this case, each mass in isolation carries a field given 
by Eq. @, e.g. 
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for the satellite mass and likewise for the parent mass 
M . For convenience, we place M at (r = 0, Z = d) and 
m at the origin (r = 0,z = 0) (see Fig. [!}:). The effective 
density becomes 



N[4>M-,4> m ] = 
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The main difference between this case and the previous 
cases is that the effective density does not have compact 
support and indeed diverges at the position of the two 
bodies. 

Nonetheless, just like for the divergent physical den- 
sities of point masses, we can still evaluate the forces 
induced by one body on the other. At the position of the 
satellite mass m 
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Since <p m is symmetric around the origin, it provides no 
acceleration there and the total is 
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FIG. 4. Numerical test of weak screening. Numerical results 
for the leading order acceleration "(0)" in the weak-screening 
appr oxim ation agree well with the analytic expression from 
Eq. {17J while those for the first order correction "(1)" from 
Eq. ( 91) indicate that the approximation is valid until the dis- 
tance to the shell s is smaller than a Vainshtein radius of the 
body r, m (vertical dashed line). These terms have opposite 
sign, implying that the acceleration ceases to grow as rapidly 
in the strong screening regime. See text for the shell param- 
eters of this example. 



such that the acceleration of m from the field of M is 



(28) reduced by 
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J/d 3 . Note that we have kept here 



the 0{r 2 ) correction to 4>a from Eq. (|4j to maintain the 
same order in r c throughout. The impact of the two-body 
interference for m « 1/ is thus muted by the fact that 
the main suppression still comes from the screening of 
the larger mass. As the two masses become comparable, 
the interference can make a more substantial reduction 
of the attractive force. 

This calculation also illustrates the role of the effective 
density in restoring momentum conservation or Newton's 
third law. By symmetry, this calculation also gives the 
acceleration of body M from the field of in 
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The sum of the two forces 
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0, (33) 



whereas without the interference the balance would be 
grossly violated for m/M <§; 1 [26] , 

IV. BEYOND WEAK SCREENING 

We expect the weak-screening regime to extend to the 
point at which the Vainshtein radii of the individual bod- 
ies is comparable to the characteristic distances in the 
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system. In our perturbative approach, the breakdown of 
weak screening is monitored by the next to leading order 
term (j)^ in Eq. (j9j). In this section we shall check to sec 
when this breakdown occurs numerically. 

For illustration purposes, we choose the shell test case 
(see Fig. [I]). Given that the effective density N vanishes 
for r > r3, this case has periodic boundary conditions 
and can be solved efficiently with fast Fourier transform 
techniques. 

As an example, we take the outer shell to be inscribed 
inside the cube of side length L = 1024 pixels: r% = L/2. 
We take the innermost shell to be of radius n = 3L/8 
and width — r\ = ri/96. We set the mass scales to be 
such that the Vamshtein scale of the shell t*m = r i/12 
and that of the mass r* TO = ri/24. Finally, we replace 
the point mass m with a constant density sphere of radius 
f m — ri/96. In Fig. |4j we show the numerical results for 
the leading order acceleration and first order correction 
given by Eq. ^ and ^ respectively. Here we plot the 
acceleration in units of the maximum Newtonian accel- 
eration for a test particle just outside of an infinitesimal 
shell of radius r\ 
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The numerical calculation of the leading order (f>^ term 



matches the analytic expression, Eq. (17) to excellent ap- 



proximation showing the discretization onto the grid does 
not cause appreciable errors. The first order correction 
<p£^ becomes comparable to the zeroth order effect when 
f*m/s ~ 1 implying that the weak-screening approxima- 
tion applies all the way to separations of a Vainshtein 
radius of m. At this point the acceleration, which would 
be zero in the absence of nonlinearity, becomes a non- 
negligible fraction of the maximal acceleration cm- 

Note that for M ^> m, the validity of the weak- 
screening approximation extends to separations much 
smaller than r*M since nonlinear corrections for M come 
in through the characteristic scale of the shell r^/Ai 
rather than the characteristic scale of the separation 
r*A//s. Thus in the shell case, the qualitative rule that 
weak screening applies until the separation becomes com- 
parable to the Vainshtein radius of the individual sources 
is too conservative if applied to both r„ n and r*M- 



With an attractive fifth force, massive objects such as 
dark matter halos in a mass shell or compensated void 
are attracted to the center of the void with an accelera- 
tion proportional to their mass. A cosmological void is 
also naturally in the weak-screening regime for the cos- 
mologically interesting case where Hqt c / '(3 ~ 1. Since 
the Vainshtein radius of a halo is comparable to its virial 
radius, the maximal effect will be on halos that are close 
to a virial radius of the edge. This effect has potential 
observable consequences for mass segregation near the 
edge of voids, but it remains to be seen in cosmological 
simulations whether deviations from the idealizations of 
spherical symmetry, perfect compensation and constant 
underdensity mask this effect. 

For the infall problem, the two-body interference pre- 
dicts a reduction of major mass mergers where the bodies 
are of comparable mass. The infall case also shows how 
the interference restores Newton's third law or momen- 
tum conservation in the joint system. 

These effects are signatures of the Vainshtein mecha- 
nism, which is itself common to massive gravity, galileon 
and braneworld scenarios. Our analytic calculations 
serve as simple illustrations that expose new aspects of 
the mechanism though constructing realistic cosmologi- 
cal tests of it will require going beyond the idealizations 
considered here. 
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V. DISCUSSION 



The weak-screening regime, where the Vainshtein 
mechanism is just beginning to suppress fifth forces, ex- 
hibits apparent or macroscopic violation of the equiva- 
lence principle. We have illustrated this effect through 
analytic calculations of the acceleration of particles 
within a mass shell, compensated void, and toward a 
spherically symmetric mass. Numerical tests show that 
weak-screening approximation is valid until the bodies 
are separated by less than a Vainshtein radius. 
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